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Cosmological implications of a stellar initial mass function that 
varies with the Jeans mass in galaxies 
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ABSTRACT 



Observations of star- forming galaxies at high-z have suggested discrepancies in the in- 
ferred star formation rates (SFRs) either between data and models, or between complementary 
measures of the SFR. These putative discrepancies could all be alleviated if the stellar initial 
mass function (IMF) is systematically weighted toward more high-mass star formation in 
rapidly star-forming galaxies. Here, we explore how the IMF might vary under the central as- 
sumption that the turnover mass in the IMF, Mc, scales with the Jeans mass in giant molecular 
clouds (GMCs), A/j. We employ hydrodynamic simulations of galaxies coupled with radia- 
tive transfer models to predict how the typical GMC Jeans mass, and hence the IMF, varies 
with galaxy properties. We then study the impact of such an IMF on the star formation law, the 
SFR- A/* relation, sub-millimetre galaxies (SMGs), and the cosmic star formation rate density. 
Our main results are: The H2 mass-weighted Jeans mass in a galaxy scales well with the SFR 
when the SFR is greater a few AIq yr^^. Stellar population synthesis modeling shows that 
this results in a nonlinear relation between SFR and Lboi, such that SFR oc Lboi°'^^- Using 
this model relation, the inferred SFR of local ultraluminous infrared galaxies decreases by a 
factor ~ x2, and that of high-z SMGs decreases by ^ x3 — 5. At z ~ 2, this results in a 
lowered normalisation of the SFR-Af„ relation in better agreement with models, a reduced 
discrepancy between the observed cosmic SFR density and stellar mass density evolution, 
and SMG SFRs that are easier to accommodate in current hierarchical structure formation 
models. It further results in a Kennicutt-Schmidt (KS) star formation law with slope of ^ 1.6 
when utilising a physically motivated form for the CO-H2 conversion factor that varies with 
galaxy physical property. While each of the discrepancies considered here could be alleviated 
without appealing to a varying IMF, the modest variation implied by assuming AfcOcMj is a 
plausible solution that simultaneously addresses numerous thorny issues regarding the SFRs 
of high-z galaxies. 

Key words: stars:luminosity function, mass function - stars: formation - galaxies: formation 
-galaxies: high-redshift - galaxies: ISM - galaxies: starburst - cosmology:theory 



1 INTRODUCTION 

The buildup of stellar mass over cosmic time is a central issue in 
understanding the formation and evolution of galaxies. A common 
approach to quantifying stellar growth is to measure the evolution 
of the star formation rates (SFRs) of galaxies. This is done using a 
wide variety of tracers from the ultraviolet (UV) to the radio. Gen- 
erally, all such measures trace the formation rate of higher-mass 
(typically O and B) stars, while the bulk of the stellar mass form- 
ing in lower-mass stars is not directly detected. Hence measuring 
the true rate of stellar growth requires assuming a conversion be- 
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tween the particular tracer flu x and the total stellar mass being gen- 
erated (e. g. iKennicuti 1 998al : ?). This requires assuming some stel- 
lar initial mass function (IMF), namely the number of stars being 
formed as a function of mass. 

On global cosmological scales, multi-wavelength observa- 
tions of galaxies are converging on a broad scen ario for the cosmic 
star formation rate evolution tMadau et al.ll996l) . Galaxies at high- 
redshift appear to be more gas-rich and forming stars more rapidly 
at a given stellar mass than present day galaxies (see the recent 
review by Shapley 2011). The cosmic star formation rate density 
rises slowly from early epochs to peak between redshifts 2: aal-3, 
and then declines toward z=0 (e.g. iHopkins & BeacomllOOq ). This 
global evolution is also reflected in measurements of the star for- 
mation rates of individual galaxies, which grow significantly at a 
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given stellar mass from today out to z ^ 2, prior to which they 
show a much slower evolution (e.g. Daddi et al.ll2005c IPavguOOSl : 



iGonzalez et al.l |201(]l : lHopkins et aljboioh . 

Meanwhile, recent advances in near-infrared capabilities have 
enabled measurements of the buildup of stellar mass out to high 
redshifts, as traced by more long-lived stars (typically red giants). 
Nominally, the integral of the cosmic star formation rate, when cor- 
rected for stellar evolution processes, should yield the present-day 
cosmic stellar mass. Analogously, the time differential of the stel- 
lar mass evolution should be the same as measured star formation 
rate. Thus in principle there is now a cross-check on the rate of 
high-mass stars forming relative to lower-mass stars. 

Preliminary comparisons al ong these lines h ave yielded gen- 
eral agreement out to z ~ 1 (e.g. lBell et alJuOOTT) . However, mov- 
ing to higher redshifts into the peak epoch of cosmic star formation, 
there are growing hints of discrepancies: the integrated cosmic star 
formation rate density seems to exceed the o bserved stellar mass 
funct i on (accounting for stellar mass loss ; Hopkins & BeacomI 



200d lElsner et al.l2008l; IWilkins et all2008l : iPerez-Gonzalez et al.1 



20081) . These discrepancies are relatively mild, at the factor of 
two to three level, so could perhaps be resolved by a more care- 
ful consideration of systematic uncertainties in SFR and A/* mea- 
sures. Indeed, som e analyses fail to show strong discrepancies (e.g. 
ISobral et al.ll2012h . Nevertheless, it is interesting that when dis- 
crepancies are seen, they unanimously favor the idea that the ob- 
served stellar mass growth appears slower than expected from the 
observed SFR. 

Theoretical models have made a number of advances to- 
ward understanding the observed properties of star- forming galax- 
ies at high redshifts. Simulations advocate a picture in which 
continual gas accretion from th e intergalactic medium (IGM) 
feeds galaxies with fresh fuel (e.g. Mo et al.l 19981: Robertson et al 



20041 : TGovemato et al.ll2009l : ICeverino et al.ll2010al: lAgertz et al 



and ultimately dr ives star for mation dKeres et al.l 120051 : 



Dekel et a l. 2009a b; Ker es et alJl200Sl) . In such models, stars typi- 



cally form at rates proportional to the global gas accretion rate from 
the IGM, regulated by feedback processes JSpringel & Hemguisli 
I2OO3I : lOave et al.l uOllh . Such a scenario predicts a fairly tight 
and nearly linear relation between star fo rmation rate a nd stellar 
mass for star-forming galaxies (e.g. lOave et al. 2000 : Finlator et al.l 
|200g). Observations ha ve now identified and quantified such a re- 
lation out to z~6-% (lNoeskeetaIll2007j |bl: lDaddietalJl2007l : 
IStark et al.l2009l : lMcLure et al.l2oTlh . which has come to be called 
the "main sequence" of galaxy formation. The bulk of cosmic star 
formation appears to occur in galaxies along this main sequence, 
while merger-driven starbursts contribut e < 20% globally (e.g. 
iRodighiero eTaLll201 ll : IWuvts et al.ll201ll) . 

Although this scenario broadly agrees with observations, a 
small but persistent discrepancy has been recently highlighted be- 
tween the evolution of the main sequence in simulations and obser- 
vations, particularly during the peak epoch of cosmic star formation 
(1 ^ -^ i$ 3). In it, the rate of stellar mass growth at these red- 
shifts in models is typically smaller by ~ 2 — 3 x than the observed 
star formation rates. This di screpancy ex ists for both cosmological 
hydrodynamic simulations (|Dava|2008l). as well as semi-analytic 
models (SAMs: lDaddi et alj2007l : lElbaz et alj2007h . and is t o first 
order independent of model assumptions about feedback dDavg 
[2008). T his is further seen in both the global cosmic star forma- 
tion rate jWilkins et al.ll2008p as well as when comparing individ- 
ual galaxies at a given M« jDavg|2008r) . Similarly, galaxy forma- 
tion simulations that utilise a variety of methods all show a paucity 
of galaxies that form stars as rapidly as galaxies with the highest 



star formation rate s at z ~ 2, the sub-millimetre galaxies ( SMGs) 
JBaugh et alJboOSi : lOave e"tani2010l : lllavward et al.||201 ih . In all 
cases, the models tend to favor lower true SFRs than implied by us- 
ing available tracers and using conversion factors based on a canon- 
ical IMF. 

One possible but speculative solution to all these discrep- 
ancies is that the stellar IMF in galaxies at z ~ 2 is different 
than what is measured directly in the Galaxy (e.g. ?). The dis- 
crepancies described above, between the various observations as 
well as between models and data, would all be mitigated by an 
IMF that forms somewhat more high-mass stars than low-mass 
ones at those epochs compared to the present-day IMIq Never- 
theless, it is important to point out that at present, no firm ev- 
idence that the IMF va ries strongly from the locally observed 
one (see the review by iBastian et al.l uOlOl) . Locally, some ob- 
servations suggest that a top-heavy/bottom-light IMF may apply 
to the Galactic C entre JNa vakshin & Sunvaev 2005; Stolteetal 
I2OO5). Similariy, iRiekeeTalJ ( 1 19931) and iForster Schreiber et al 
( 12003,) suggest a tu rnover mass a factor of ~ 2 — 6 larger 
than in a traditional iKroupal j2002t) IMF in the nearby starburst 
galaxy M82. Simultaneous fits to the observed cosmic star forma- 
tion rate density, integrated stellar mass measurements, and cos- 
mic background radiation favour a "paunchy" IMF that produces 
more stars at intermediate masses (IFardal et al.l2007l) . lvan DokkumI 
(2008) suggested that the IMF may be more top-heavy at high- 
redshift {z ~0.8) based on an analysis of the evolution of the 
colours and mass-to-light ratios of early-type galaxies. However, 
these observati ons can all be inte rpreted without the need for 
IMF variations ( IBastian et al.ll201CI) . Beyond this, some observa- 
tions fi nd evidence for a bottom-heavy IMF in z=Q earl y-type 
galaxies('van Dokkum & Conrovl2oTTI : ICappellari et al.l2012 : ?). It 
is therefore interesting to examine whether an IMF-based solution 
is viable and consistent with a broad suite of observations, both 
locally and in the distant Universe. 

In this paper, we explore the cosmological consequences of a 
physically-based model for IMF variations. Past work has generally 
focused on empirically determining the amount of IMF variation 
needed in order to solve one (or more) o f the above problems (e.g. 
Fardal et al.l l2007l : Ivan DokkumI I2OO8I; lDavdl2008l : Iwilkins et al.) 
2008). Here, instead, we make a single critic al assumption, first 
forwarded by leans, and later expanded upon bv LarsonI (120051) and 
?: The IMF critical mass (Mc) scales with the Jeans mass in a gi- 
ant molecular cloud (GMC). For reference, we call this the leans 
mass conjecture. We employ hydrodynamic simulations of isolated 
galaxies and mergers including a fully radiative model for the inter- 
stellar medium (ISM) to predict the typical Jeans mass of GMCs in 
galaxies with different physical conditions, corresponding to qui- 
escent and starbursting systems both today and at z = 2. Applying 
the Jeans mass conjecture, we then make a prediction for how the 
IMF varies with global galaxy properties, and explore the implica- 
tions for such variations on the discrepancies noted above. 

We emphasise that the main purpose of this paper is to utilise 
numerical models of the molecular ISM in galaxies to investigate 
the consequences of an IMF in which Mc scales with A/j . We do 
not directly argue for such a scaling relationship; this is taken as an 
assumption, and we only seek to study its implications in a cosmo- 



^ This could be described as a "top-heavy" IMF, which we specify as an 
IMF whose high-mass slope is different than local, or a "bottom-light" IMF, 
which we define as retaining the same high-mass slope but forming fewer 
low-mass stars. This paper focuses on bottom-light IMFs. 



Implications of an IMF that Varies with the Jeans Mass in Galaxies 3 



logical context. For the reader's edification, we present arguments 
both for and against the Jeans mass conjecture in §[5] 

This paper is outlined as follows. In §[2l we detail our numer- 
ical models. In §[3] we discuss variations of the IMF versus galaxy 
star formation rates. In § [4] we investigate the effect of this model 
on derived star formation rates, paying particular attention to the 
star formation law (ij 14. H and §|4J2f, the SFR-M* relation at high- 
redshift (i; 14.3b . and the evolution of the cosmic star formation rate 
density (§|43}. In §[5] we present a discussion, and in §(6] we sum- 
marise. 



2 NUMERICAL METHODS 

Our main goal is to simulate the global physical properties of 
GMCs in galaxies, and understand the effect of the varying 
Jeans mass on observed star formation rates. These methods and 
the corresponding equat ions are described in significant detail in 
[Narayanan et al.l ( l2011bl lcl): for the sake of brevity, we summarise 
the relevant aspects of this model here, and refer the interested 
reader to those papers for further detail. 

We first simulate the evolution of galaxies hydrodynami- 
cally using the publi cly available code GADGET- 3 JSpringefcoOSl : 
ISpringel et al.ll2005n . In order to investigate a variety of physical 
envirormients, we consider the evolution of disc galaxies in iso- 
lation, major (1:1 and 1:3) mergers, and minor (1:10) mergers at 
both low (z=0) and high (z =2) redshifts. These sim ulations are 
sunmiarised in Table Al of iNaravanan et al.l ( 201 id) . The discs 
are initialised according to the IMo et al.l rl9981 model, and em- 
bedded in a live dark matter halo with a[Hemquist ( 1990) density 
profile. Galaxy mergers are simply mergers of these discs. The halo 
concentration and virial radius for a halo of a given mass is moti- 
vated by cosmological Ai'-body simulati ons, and scaled to m atch 
the expected redshift-evolution following iBuUock et al J ( 120011) and 
[Robertson et al. (2006). 

For the purposes of the hydrodynamic calculations, the ISM 
is mo deled as multi-phase, with cold clouds embedded in a hotter 
phase dMcKee & OstrikeilI977h . The phases exchange mass via ra- 
diative cooling of the hot phase, and supernova heating of cold gas. 
Stars form within the cold gas according to a volumetric Kennicutt- 
Schmidt star formation relation, SFR ~ Prn M. with normalisa- 
tion set to ma t ch th e locally observed relation ( lKennicuttlll998al : 
ISDringelll2000l: ICoxet al. 2006). 

In order to model the physical properties of GMCs within 
these model galaxies, we perform additional calculations on the 
SPH models in post-processing. We first project the physical prop- 
erties of the galaxies onto an adaptive mesh with a 5'^ base, span- 
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ning 200 kpc. The cells recursively refine in an oct-subdivision 
based on the criteria that the relative density variations of metals 
should be less than 0.1, and the V-band optical depth across the 
cell TV < 1. The smallest cells in this grid are of order ~ 70 pc 
across, just resolving massive GMCs. 

The GMCs are modeled as spherical and isothermal. The 
H2-HI balance in these cells is calculated by balancing the pho- 
todissociation rates of H2 aga inst the growth r ate on dust g rains 
following the rnethodol ogy of iKrumholz et alj ( |2008| . l2009ah and 
iKrumholz et alj (2009b). This assumes equilibrium chemistry for 
the H2. The GMCs within cells are assumed to be of constant 
density, and have a minimum surface density of 100 Mq pc~^. 
This value is mot ivated by the typical su r face density of lo- 
cal group GMCs dMcKee & Ostrikej l2007l : iBolatto et al.] l2008l : 
iFukui & Kawamurall20IOt) . and exists to prevent unphysical condi- 
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Figure 1. Evolution of SFR, enhancement in Af j (above 250 Mq , a typical 
Jeans mass for an n = lO'^cm"'^, T = 10 K GMC), H2 mass-weighted 
GMC temperature and H2 mass-weighted GMC density in a model 2=0 
disc galaxy, 2=2 disc galaxy, and 2=0 1:1 major galaxy merger. During the 
merger, warm dust heated by the starburst couples efficiently with the gas 
and raises the temperature of H2 molecules. This drives a rise in the typical 
Jeans mass of GMCs. 2=0 quiescent discs, on the other hand, have much 
more moderate conditions in their ISM. High-2 discs represent intermediate 



tions in large cells toward the outer regions of the adaptive mesh. In 
practice, the bulk of the GMCs in galaxies of interest for this study 
(i.e., galaxies where the galactic environment has a significant im- 
pact on the Jeans mass) have surface densities above this fiducial 
threshold value. With the surface density of the GMC known, the 
rad ius (and consequently me an density) is known as well. Follow- 
ing |Narayanan_et_alJ(l20ilb|), to account for the turbulent compres- 
sion of gas, we scale the volumetric densities of the GMCs by 
a factor e'^p^^, where ap « ln(l + Mi^n /4), and Mip is the 
1 dimensional Mach num ber of turbulence dOstriker et al.luOOll : 
IPadoan & Nordlundll2002l) . 

Because this study centres around understanding the Jeans 
scale in GMCs, the t hermal state of the mo le cular ISM must 
be kno wn. Following Krumholz et alj ( 1201 Ibf) : [Narayanan et al.l 
( 1201 Ibl) and lNaravanan et al.l ( 120110) . we consider the temperature 
of the H2 gas as a balance between heating by cosmic rays, grain 
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photoelectric effect, cooling by CO or CII line emission, and en- 
ergy exchange with dust. The cosmic ray flux is assumed to be that 
of the mean Galactic valuqj. Formally, if we denote heating pro- 
cesses by r, cooling by A, and energy exchange with ^, we solve 
the following equations: 



r. 



+ TCR - Aline + *gd 
Tdust — Adust — *gd ■ 







(1) 

(2) 



We refe r the reader to Krumholz et al.ldlOl Ibh and lNaravanan et alj 
( l2011b» for the equations regarding the photoelectric effect and 
cosmic ray heating terms. 

The dust temperature is calculat ed v ia the pub licly available 
dust radiative transfer code SUN RISE djon sson 2006: |jonsson et al.l 
l201u Ijonsson & Prim ack 2010). We consider the the transfer of 
both stellar light from clusters of stars as well as radiation from 
a central active galactic nucleus (though this plays relatively lit- 
tle role i n these sim ulations). The stars emit a STARBU RST99 
spectrum JLeitherer e t al. 1999; IVazquez & Leithereri2005l) . where 
the metallicity and ages of the stellar clusters are taken from 
the hydrodynamic simulations. The radiation then traverses the 
galaxy, being absorbed, scattered and reemitted as it escapes. 
The evolving dust mass is set by assuming a constant dust to 
metals ratio comparable to the mean Milky Way value JDwekl 
ri998l : IVla dilo 1998; Calura et al. 2008), and takes the fom of the 
R = 3.15 .Weingartner & Draine. C2001.) grain model as updated 
bv lOraine & Lil ( 120071) . The dust and radiation field are assumed to 
be in equilibrium, and the dust temperatures are calculated itera- 
tively. Energy exchange between gas and dust becomes important 
typically around n ~ lO^cm"^. Henceforth, we shall refer to this 
as the density where grain-gas coupling becomes important. 

The gas cools either via CII or CO line emission. The frac- 
tion of H2 where most of the carbon is in the form of CO (versus 
atomic form) is det eiTnined following the semi-analytic model of 
IWolfire et al.l ( 1201(1) , which is a metallicity dependent model (such 
that at lower metallicities, most of the carbon is in atomic form, and 
at higher metallicities, it is mostly in the form of CO): 



-4(0. 53-0. 0451n 

/CO = /h2 X e "h/c 



-0.0971nZ')Mv 



(3) 



where Go is the UV radiation field with respect to that of the 
solar neighbourhood, and Ay is the visual extinction. Ay ~ 
i.87xi'^"cm--.^ ^' ' IWatsonl l201ll). If /co > 0.5, we assume that 
the gas cools via CO line emission; else, CII. The line emission 
is calculat ed via the publicly available es cape probability code as 
detailed in lKrumhoIz & Thompson! (20071). 

While the temperature model is indeed somewhat compli- 
cated, the typical temperature of a cloud can be thought of in terms 
of the dominant heating effects at different densities. At low densi- 
ties (n ~ 10 — 100 cm~^), the gas cools via line cooling to ~ 8 
K, the characteristic temperature imposed by cosmic ray heating. 
At high densities (n > 10* cm^'^), grain-gas coupling becomes 
efficient, and the gas rises to the dust temperature. At intermediate 
densities, the temperature is typically in between 8 K and the dust 
temperature. 



^ Tests described by iNaravanan et al.l J2011ch show that even under the 
assumption of a scahng where cosmic ray flux scales with the star formation 
rate of a galaxy, the typical thermal profile of H2 gas in a given galaxy is 
unchanged. This is because, in these environments, energy exchange with 
dust tends to dominate the temperature. 



3 THE RELATION BETWEEN JEANS MASS AND STAR 
FORMATION RATE 

We begin by examining the conditions under which the Jeans mass 
in GMCs in a galaxy may vary. We define the characteristic Jeans 
mass, 7l/j, as the H2 mass- weighted Jeans mass across all GMCs in 
a galaxy, and consider the deviations of this quantity. The density 
term in the Jeans mass calculation is the mean density of the GMC. 
In Figure[T] the top panel shows the SFR as a function of time for an 
unperturbed z=0 disc galaxy, a z=0 major 1:1 galaxy merger, and a 
z=2 disc. The second panel shows the evolution of A/j/250 Mq. 
250 Mq is the Jeans mass for physical conditions as found in local 
discs like the Milky Way (MW), namely n = 1 — 2 x lO^cm"^ 
and T = 10 K. Our disc galaxy has a baryonic mass of ^ 5 x 10^" 
Mq, and is the fiducial MW model in Nara yanan et alj ( l201 Ibl) . We 
note that the fiducial z=0 MW model naturally produces densities, 
temperatures, and Mj as observed in the MW, thereby providing an 
important check for our ISM model. The major merger is a merger 
of two A/bar = 1-5 X 10^^ Mq galaxies. The high-z disc model 
is an unperturbed Afbar = IO^^Mq disc. The remaining panels 
show the evolution of the characteristic temperatures and the char- 
acteristic GMC densitieo 

For the bulk of the z=0 disc's life, and during non-interacting 
stages of the galaxy merger, the GMCs are relatively quiescent, re- 
taining surface densities near 100 Mq pc^'^, and temperatures near 
the floor established by cosmic ray heating of ~ 10 K. During the 
merger-induced starburst, however, the story is notably different. 
Gas compressions caused by the nuclear inflow of gas during the 
merger cause the average density to rise above lO^'cm^'^ in GMCs. 
At this density, grain-gas energy exchange becomes quite efficient. 
At the same time, dust heating by the massive stars formed in the 
starburst causes the mass-weighted dust temperature to rise from 
~ 30 K to ~ 80 K. Consequently, the kinetic temperature of the 
gas reaches similar values. While the mean gas density also rises 
during the merger, the Jeans mass goes as r^"/n^". As a re- 
sult, Mj rises by a factor ~ 5 during the merger-induced starburst. 
The z=2 disc represents an intermediate case: high gas fractions 
and densities drive large SFRs, and thus warmer conditions in the 
molecular ISM. 

The characteristic Jeans mass in a galaxy is well- 
parameterised by the SFR. In Figure |2] we plot the star formation 
rates of all of the galaxies in our simulation sample against their 
A/j. The points are individual time snapshots from the different 
galaxy evolution simulations. At high SFRs, dust is warmed by in- 
creased radiation field, and thermally couples with the gas in dense 
regions. The increased temperature drives an increase in the Jeans 
mass, resulting in roughly a power-law increase of Afj with SFR. 

This relation does not extend to SFR J; 3 — 5 MQyr~^. 
Here, the mean densities become low enough (n < lO^cm^'') that 
energy exchange with dust no longer keeps the gas warm. Cos- 
mic rays dominate the heating in this regime, and the gas cools to 
~ 8 — 10 K, the typical temperature that results from the balance 
of CO line cooling and cosmic ray heating. Because the tempera- 
tures are relatively constant, and the characteristic densities of or- 
der ~ 10 — 100 cm~^, the Jeans mass flattens with respect to SFR. 
There may be additional variations in Af j at smaller masses, owing 
e.g. to metallicity variations that we neglect here. For our puiposes, 
we are predominantly interested in higher- 2 galaxies that generally 



•^ Note that because Mj is calculated as the H2 mass-weighted Jeans mass, 
one cannot convert from the characteristic temperatures and densities shown 
in Figure[T]to Afj. 
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have SFR^SMQyr '^ , so we do not try to characterise in detail the 
behavior of Mj at lower SFR. 

The resulting relation is well-fit by a powerlaw 
Mj/250Mq = 0.5 X SFR"'' at SFR > SMQyr'^ This 
fit is shown in Figure |2] by the black line, and extends up to 
'-^ 1000 Mq /yr in our models. The fit excludes the low SFR 
tail at SFR < 3 MQyr^^. There is substantial scatter around the 
relation owing to the fact that at a given SFR, there may be a 
variety of physical conditions. For example, galaxies with SFR 
~ 50 MQyr~^ may represent merger-induced starbursts with 
extremely warm nuclear star formation, as well as high- 2 discs 
with more distributed star formation and moderate temperatures. 

We choose to parameterise A/j in terms of the SFR in order to 
facilitate analysis of inferred observed SFRs of high-z galaxies. In 
principle, one could imagine a relationship between Mj and stellar 
mass. Indeed, broadly, the lower stellar mass galaxies tend to have 
low Af J, and vice versa. However, due to similar stellar masses on 
(e.g.) before and after a starburst event, the scatter in an Afj-Af» 
relation is even larger. Similarly , AJfj may correlate well with the 
H2 gas surface density. Indeed, [Narayanan et al.l ( 1201 Iq) showed 
that the physical conditions in the gas (temperature and velocity 
dispersion) scale well with the molecular gas surface density. The 
scatter in such a relation, however, is similar to the Afj-SFR rela- 
tion (the SFR is forced to scale with the volumetric gas density, and 
shows a tight correspondence with th e gas surface density for these 
simulations: lNaravanan et al.l2011an . 

One could similarly imagine a Jeans mass that is determined 
by the median density in the gas, rather than the mean density. As 
it turns out, this has little effect on our results. We discuss this issue 
in further detail in Appendix [A] 

Finally, we emphasise that what governs the Jeans mass in 
GMCs is the physical conditions in the ISM, not the global mor- 
phology. As such, both heavily star-forming discs (e.g. the BzK 
population at z ~ 2) as well as galaxy mergers may exhibit warmer 
ISM gas relative to today's discs. One can expect short-lived nu- 
clear starbursts driven by major mergers in the local Universe to 
exhibit an increase in Afj by a factor of > 5. Similarly, high-z 
disc galaxies that are forming stars at many tens of MQyr"'^ for a 
much longer duty cycle may show an increased Afj in their GMCs 
by a factor '^ 2 — 3. The increase in temperature of the ISM is 
typically driven by warm dust and increasing fractions of molecu- 
lar gas above the gas-grain coupling density with increasing SFR. 
We note that very high redshifts, when comparing to galaxies at 
a fixed density, then the Jeans mass will only scale with the SFR 
once the mean gas temperature is warmer than CMB temperature. 
We discuss this more quantitatively in i) |4.5l 



4 IMPLICATIONS FOR OBSERVED STAR FORMATION 
RATES 

We showed that the H2 mass-weighted Jeans mass can vary sub- 
stantially in galaxies with physical conditions different than that in 
the Milky Way. We now employ make the assumption that varia- 
tions in A/j directly into variations in the characteristic mass A/c 
of the IMF. This then allows us to investigate how such variations 
in Mj will change the inferred star formation rates. 

It is useful to parameterise this variation in terms of the how 
the bolometric luminosity Lboi varies with SFR. In a star-forming 
galaxy, Lboi is essentially dominated by the emission from high 
mass stars. All currently used extragalactic SFR tracers trace such 
stars. Hence examining the variations in Lboi is likely to be broadly 



^ 



000 


■ • z=0 Discs 
■ High-z Discs 

. A z=0 Mergers 
+ High-z Mergers 




#<^^^' 


100 





10 



100 

SFR (Meyr ') 



1000 



Figure 2. Relationship between SFR and H2-mass weighted Jeans mass of 
GMCs for all simulation snapshots in our sample. The points show individ- 
ual simulation snapshots, and are colour-coded by the type of simulation. 
The solid line denotes the best fit power law, while the dashed line shows 
a fiducial value of 250 Mq, the typical Jeans mass for a GMC with den- 
sity n fsj 1 — 2 X 10^ cm^'^, and temperature T Si 10 K. Because dust 
and gas heating is dominated by star formation, A/j scales with the SFR. 
At low SFRs, the mean GMC density drops enough that gas-grain coupling 
becomes inefficient, and the heating is dominated by cosmic rays. In this 
regime, the temperatures of the GMCs are roughly constant (~ 10 K), as 
ai'e the Jeans masses. 



representative of the variations that would be seen in any of the 
canonical SFR tracers. For an invariant IMF, Lboi is linearly related 
to SFR (assuming, as we will throughout, that there is no significant 
contribution from active galactic nuclei). But when the IMF varies, 
then Lboi will trace only high-mass stars, while the true SFR in- 
cludes all stars, thereby breaking the linear relation. 

To quantify the SFR-Lboi relation, we begin by assuming that 
the IMF has a broken power-law form characterised by a turnover 
mass, A/c, a low-mass slope of 0.3, and high-mass slope of -1.3 
(where the IMF slope specifies the slope of the logio(dN/dlogM)- 
logio(M) relation). If Afc scales with A/j, then we can calculate 
the variation in the Lboi-SFR relationship as a function of Mc by 
employing stellar population synthesis modeling. 

We e mploy FSPS, a publicly available stellar populatio n syn- 
thesis code JConrov et al.l2009ll20I0l : IConrov & Gunnl2010l) . fsps 
allows for the quick generation of stellar spectra given a flexible in- 
put form for the IMF. We assume a constant star formation history 
(SFH) of 10 MQyr"^. This was chosen as a reasonable approxi- 
mation for the star formation history ofa_(z=0)Af* — IO'^'^'^Mq 
galaxy in cosmological simulations ( lDavg|2008h . We additionally 
assume solar metallicity stars, and include the contribution from 
thermally pulsating asymptotic giant branch (AGB) stars. No dust 
extinction is included. 

Using FSPS, we find that Lboi/SFR 
model assumptions, we vary the assumed star formation history 
and assume a constant SFH of 50 MQyr~^, as well as an exponen- 
tially declining SFH with an initial SFR of 50 and 10 MQyr^^ and 
an e-folding time of 10 Gyr. All assumed SFH models return an 
Lboi/SFR-A/c exponent within 10% of our fiducial model. Com- 
bining our fiducial Lboi/SFR-Afc relation with the model fit to the 
SFR-Afc relation fit from Figure|2] we arrive at the relation: 
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Figure 3. Kennicutt-Schmidt Ssfr 

line shows relation Ssfr 

ing the SFR corrections derived for the model IMF in this paper, and 

th e model function a l form for the CO-H2 conversion factor presented 

in [Narayanan et alj bOlla) . The blue dashed line denotes the relation 

Ssfr ~ ^moP which results if the IMF does not vary with the Jeans mass 

of GM Cs. The low-z data are from iKennicutt et al.1 00071); [Wong & Bhtz 



E 



2002 ); JCrosthwaite & Turned ( l2007b ; ISchuster et al.l ( l2007b. Ifiigiel et al 

200s ); lKennicutli V1998b[). whereas the high-z data are from faenzel et al 

2O10l):lDaddi et alH20I0j|bl); Bothwell et al. 1 2010); Bouche et alj to07l) : 

et al.l J2005l) ; lTacco'niet alj l l200ai2008l) ; lEngel et al.ll2O10r . 



SFR = 



Lb 



10iOLf7 



'0 



Mgyr^^ , for SFR > SMgyr^ 



(4) 



As seen in Figure |2] at SFR < 3MQyr ^ tiiere is no strong varia- 
tion of Mj with SFR, so in that regime SFRcx Lboi- 

At ULIRG-type SFRs (~ 100 MQyr"^), our model sug- 
gests that SFRs may need to be adjusted by a factor ~ 2. 
At the luminosities of distant SMGs, the most luminous, heav- 
ily star-forming galaxies in the Universe, typical SFRs may be 
closer to ~ 400 MQyr~^, rather than the canonically inferred 
~ lOOOM0yr"\ 

Though the correction to the inferred star formation rates from 
a standard Kroupa IMF is relatively mild (only a factor of ~ 3 at 
a SFR of ~ 1000 MQyr~^), the impact on cosmological star for- 
mation rates can be substantial owing to the strong evolution of 
global SFRs. In the following subsections, we discuss the impact 
of a SFR-Lboi relation that derives from a variable IMF on the ob- 
served Kennicutt-Schmidt star formation relation, the SFR-i\f, re- 
lation, high-redshift sub-millimetre galaxies, and the observed cos- 
mic star formation rate density. 

4.1 The Kennicutt-Schmidt Star Formation Law 

The Kennicutt-Schmidt (KS) star formation relation is a power- 
law relationship between the star formation rate surface density 
against the gas surface density of both sub-regions of galaxies, 
as well as galaxies as a whole. A reasonably tight re lationship 
has been observed on scales as small as GMC clum ps jWu et al.l 
l2005ll2010l;lHeiderman et alJlOldlLada et al.l201 ll), ~ 10^ - 10^ 



pc sc a le regions in galaxie s (e.g. iKennicutll 



1998bl : iBigiel et al.1 



low and high-z (e.g. JOao & Solomon! l2004al|bl; Narayanan et al.l 
12005; Bussman n^ alj|2008l; llono et alj|2009l;paddi et alj|2010ij ; 
loenzel et alj I20ld ; iGarcia-Burillo et al.ll201l() . The relationship 
between the star formation rate and gas density is crucial both 
for our understanding of the dom inant physi c al effects in con- 
trolling star formation ra t es (e.g. ISilkl Il997t lElmegreenI 120021 ; 
iKrumholz & McKee|[20o3 ; lTaiJl201Cl ; lOstriker & Shett^^l201 ll. to 



name just a f ew), as well as for inp ut into galaxy e volution sim- 
ulations (e.g. [ Springel & Hemguistl 12003; Robertson & Kravtsovl 



l2008l ; lBoumaud et alj201(]|;ICeverino et al.l2010bh . 



The KS relation is typically measured observationally via a 
combination of inferred star formation rates, and measured H2 
gas masses. At face value, assuming a constant Lboi-SFR rela- 
tion, and a constant CO-H2 conversion factor, a relation Esfr — 
E^Vo"^' results when including both local and high-z data (e.g. 

However, it is now well 



^PT, iGauii.a vym^ii iiii^luding both lo 

iDaddi et allboiOal ; iGenzel et alJl201oH^ 



established tha t the CO-H2 conversion factor is not universal. Ob- 
serva t ions (e.g. Bolatto et al.l2008l ; lTac coni et al. 2008; Ler ov et al.l 
201ll;lGenzel et al.ll201lh. as well as theory (.Shettv et al...2011b 



Feldmann et al.l201 IclNaravanan et al.l201 Ibl ld) motivate a contin- 
uous form of Xco in ter ms of the physical properties of the H2 gas. 
[Narayanan et al.l ( 120110) derived a functional form for the CO-H2 
conversion factor (Xco) as a function of CO surface brightness 
as an observable proxy for the gas density and velocity dispersion. 
These authors showed that combining their model of Xco with lit- 
erature data results in a steeper inferred KS relation, Esfr ^ ^^2- 
This owes to decreasing CO-H2 conversion factors at higher SFR 
surface densities, and assumes a linear Lboi-SFR relation. 

If the stellar IMF varies with the physical conditions in the 
ISM, the KS relation will be further modified, since Lboi-SFR is no 
longer linear Equation|4]suggests when utilising a standard Kroupa 
IMF, the inferred SFR for the most heavily star-forming galaxies 
may be overestimated, and that the KS power-la w index will de- 
crease. In Figure[3] we plot the compiled data from lNarayanan et al.l 
(2011c) which includes both resolved regions from nearby galfix- 
ies, as well as global data from quiescent discs and starburst galax- 
ies from 2: « — 2, on a KS relation. The data have had their SFR's 



calculated via Equation B] The best fit relation is Ssfr 
reference, we show the Narayan an et al.l (l2011cl) relati 



relation. 



^H2-r 

Ssfr 



I2OO8I ; ISchrubaetai] bouh . and galaxies as a whole both at 



Eh2 relation via the blue dashed line. 

In order to calculate the modified SFRs in Figure|3] we assume 
that the original SFRs reported in the papers that the data origin ate 
from are relateable to the Lir via the standard JKennicuttl <1998al) re- 
lation, and that Lboi ~ Lir. This is likely a reasonable assumption. 
Galaxies below SFR of 3 MQyr^^ do not undergo a modification 
to the SFR in our model. The galaxies which will have the strongest 
modification to their SFR when utilising our model IMF are those 
that form stars at the rates of SFR > 100 MQyr^^. These galaxies 
are typically mergers in the local Universe, and either heavily star- 
forming discs or mergers at high-z, and the bulk of their luminosity 
is in the infrared. 

Any theory that relates the SFR to the free-fall time in a GMC 
naturally results in a relation SFR ~ pjj'l . For galaxies of a constant 
scale height, this translates into a similar surface density relation. 
An assumption of a constant scale height across such a wide range 
of physical conditions is debatable, however (??). Alternatively, if 



* We note that when considering only local data from qu iescent discs, the 
KS relation has a roughly hnear index JBigiel et alj|2008l) . It is only when 
including higher surface-density galaxies such as ULIRGs or high-z g alax- 
ies that the fit is torqued toward higher values ( IKrumholz et alj2009br) . 
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Figure 4. Similar to Figure |3] but with gas surface densities divided by the 
dynamical time (top) or free-fall time (bottom). Both relations are nearly 
linear, with star fomiation efficiencies eg ~ 0.01. See text for details. 



SFR ~ p for quiescent regions, and has a steeper relation for star- 
burst galaxies as ex pected for a scenario in which gas dominates 
the vertical gravity (lOstriker & ShettvlbOllb . then a surface den- 
sity relation steeper then unity is natural. 



4.2 The Dynamical Time Star Formation Law 

As an alternative to the standard KS relation that relates Esfr to 
Eh2, some theories posit that the gas density divided by the galaxy 
dynamical time or local free fall time may be the relevant physical 
parameter controlling the SFR (e.g. ISilkl 1 9971: lElmegreenl2002l) . It 
is of interest, therefore, to estimate the modified Esfr — EH2/idyn 
and Esfr — EH2/iff relations given our model IMF which varies 
with the physical conditions in the ISM. We plot these relations in 

Figure |4] 

Following iDaddi et al.l J2010ah and ICenzel et al.l ilOld) . we 
define the dynamical time of the galaxies in Figure [3] as the ro- 
tational time either at the galaxy's outermost observed radius, or 
the half light radius, depending on the sample. When assuming 



a constant SFR-Lboi relatioro Es fr is nearly linearly re lated to 
Emoi/idyn, with expoueut ~ 1.03 JNaravanan et al.l201 id) . When 
we use Equation |4] to calculate the star formation rates, a relation 
Esfr ~ Emoi/idy^n results. Hence this IMF-driven variation pro- 
duces only a minor change to the relation relative to a constant 
SFR-Lboi relation. 

Why does this relation change only marginally when the 
change in the Esfr — Emoi was more pronounced? Emoi/idyn 
spans ~ 6 orders of magnitude when considering the galaxies in 
Figure |3] whereas Emoi spans only ~ 4. The drop in the SFR of 
the most heavily star forming galaxies when assuming our model 
IMF makes less of an impact on the inferred relation than in the 
Esfr - Emoi case. 

Similarly, we can consider how our model affects the Esfr — 
Emoi/iff relation. Both observations and models have suggested 
that the star formation efficiency per free fall time should be 
roughly constant in molecular gas (Kmmholz & McKee 20Qa; 
Krumh olz & Tan" 20071: Isvans et al. 2009: Ostriker & Shett vlbOlTl : 
Krumholz et aL201 laO . Following . Krumholz et al. 1.2011a ). we in- 



fer the free fall time from the observable properties of the galaxy. 
Similar to the Esfr — Emoi/idyn relation, when assuming a 
constant SFR-Lboi relation, Esf r varies roughly linearly with 
Emoi/iff JNaravanan et al.ll201 id) . Introducing a variable IMF in 
the SFR calculations results in a moderate change the Esfr — 
Emoi/iff exponent to ~ 0.9 for similar reasons as in the tdyn cal- 
culation. A star formation efficiency eff ~ 0.01 provides a good fit 
to the relation Esfr = eff Emoi/te. 



4.3 The SFR-A/, Relation 

The stellar masses and star formation rates of galaxies have 
been observed to have a relatively tight relationship as early as 
z=6-S (Noeske et a l. 2007a, b: Daddi et al. 2007: Stark et al. 200^ 
iMcLure et al.ll201 11) . The SFR-i\f, relation appears to steadily de- 
crease in normalisation with time, in accordance with cosmologi- 
cal simulations that suggest that star formation is fueled by grav- 
itatio nal infall of baryons into the galaxy (e.g. Keres et akluOOSl : 
Dave etalj 1 2011 ). which shows a rapid drop with redshift (e.g. 
iDekel et alj2009bl) . 

As alluded to in § [T] while simulations predict the slope and 
scatter of the SFR-M* relation with reasonable success, the ampli- 
tude at 2 ~ 2 is less well reproduced. In models, the SFRs at a given 
M, at z — 2 are too low by x 2 — 3 compared with observations. 
This is a fairly small discrepancy given the systematic uncertain- 
ties, but is nonetheless persistent for simulations performed via a 
variety of techniques (hydrodynamic simulations performed using 
SPH or adaptive mesh refi nement, as w ell as SAMs). 

This discrepancy led .Davg ( 120081) to suggest that perhaps the 
SFRs of 2: ~ 2 galaxies are overestimated. Dave empirically de- 
rived an IMF that varies its critical mass with redshift as Mc = 
0.5(1 + z)^ Mq for < 2 < 2 as being able to reconcile the mod- 
els an d data. Our model for a varying IMF differs from the lOavg 
OOOSi) model in that here, the IMF varies with the temperature and 
density of GMCs in the ISM, rather than with redshift. Thus, a 
heavily star-forming galaxy in the local Universe may have an IMF 
more comparable to a luminous disc at z ~ 2 than a present-epoch 



^ Note, that this relation depends on the usage of a physically motivated 
functional form for the CO-H2 conversion factor that depend s on the 
CO su rface brightness and the metallicity of the galaxy CNaravanan et alj 

bond) . 
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Figure S. SFR-A/* relationship for 2=0,1,2 galaxies. Left columns are literature val ues of SFR, and rig ht show corrected SFRs using our model IMF. In the 
top two ro ws, the shaded grey region shows the median observed relationship from lElbaz et al.l | |2007|). with dispersion , and in the bottom row the crosses 
come from lDaddi et al.l 120071) . In all panels, the circles with error bars show the simulated points from Pave et al.l ( 120 111 ) in bins of M*, with dispersion. The 
red lines in the right column show the best fit to the modified SFR- .A/, relation. When using a standard Kroupa IMF, the simulated galaxies underpredict the 
SFRs by as much as a factor of ~ 10. When assuming an IMF that scales with the thermal conditions in the molecular ISM, the simulated SFR-A/, relation 
and observed one comes into much better agreement. This is most apparent in the z=2 case (bottom row). 



quiescent disc. Nevertheless, the impact of such a varying IMF does 
imply a strong dependence in redshift because high-redshift galax- 
ies form stars far more rapidly. 

In Figure |5] we plot the observed SFR-Af» data from the 
Sloan Digital Sky Survey (SDSS; z=Q) and Great Observatories 
Origins Deep Survey (GOODS; z=l-2), overlaid by the predicted 
SFR-Af , relation from simulated ga laxies drawn from the hydrody- 
namic cosmological simulations of Dave et alJoOlli) . The left col- 
umn shows the observed and modeled SFR-Af, relation assuming 
a Kroupa IMF with turnover mass Afc=0.5 Mq. With increasing 
redshift, the discrepancy between the normalisation of the observed 
SFR-Af, relation and modeled one grows, becoming noticeable at 
z> 1- 

If the turnover mass in the IMF varies with the Jeans mass in 
GMCs, this discrepancy may be ameliorated. Turning now to the 
right column, we show the same comparison after correcting the 
SFRs via Equation|4] While the relatively quiescent galaxies in the 
observed low-z samples undergo relatively little correction to their 
already modest SFRs, the galaxies at high-z tend to be brighter at 
a given stellar mass than their low- 2 counterparts due to increased 
gas fractions and densities. The SFRs for these galaxies undergo a 



larger coiTection per Equation |4] and bring the observed SFR-Af, 
relation into agreement with simulations. 

We fit a power-law form to observations, as reinterpreted using 
our evolving IMF. For the z= Q and z=\ da t a, we fit the change in the 
median relations reported bv Elbaz et al.l ( |2007[). For the z=2 da ta, 
we perform a polynomial fit to the data from lDaddi et al.l (|2007|jj. 
Our best fit relations are: 



(^ = 0) : SFR(M0yr-i) = 3.9 x 
(^ = 1) : SFR(MQyr-^) = 19.3 x 
(2 = 2): SFR(MQyr-i) = 50.2 x 



Af. 



1O"M0 
1O"M0 



IQiiM, 







(5) 
(6) 
(7) 



and are denoted in Figure|5]by the solid red lines in the right col- 
umn. 



® Note that we do not have error bars available for the iDaddi et alj 420071) 

data. 
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It important to note that all of these relations should be re- 
garded as tentative - they represent relatively limited ranges of 
stellar masses, and for the z=0 and z=l data the relations were 
derived by re-fitting published fits, not the actual data with error 
bars. This is therefore only intended to be illustrative of the sense 
and magnitude of the changes in SFR-M* when applying an IMF 
wi th Mc PC Mj. The best fit exponents for 2 :=[0,1,2] as reported 
bv'Paddi etai] J200 7') and'El baz et al] ilOOlh are ~ [0.77,0.9,0.9] 
respectively, which in our scenario decrease to ~[0. 65,0.77,0. 55]. 
This of course owes to the fact that in our model, galaxies with 
larger inferred SFRs (using a standard Salpeter or Kroupa IMF) 
undergo a larger correction to their SFR when the IMF varies with 
the Jeans mass. Similarly, the amplitudes of the SFR- A/, relations 
in Figure|5]are lower than previously estimated due to lower SFRs. 

The cosmological simulation predictions now match the ob- 
servations quite well in amplitude (around ~ L*), but the slope 
is now more discrepant. A selection effect that could explain this 
discrepancy is that low SFR galaxies tend to be absent from these 
samples, creating a floor in observed SFRs that increases to higher 
redshifts. Combined with the (not insubstantial) scatter in the SFR- 
M« relation, this could induce a shallower observed slope relative 
to the true slope (N. Reddy et al., in prep.). On the other hand, it 
could also just indicate a failing in these simulations. We leave a 
more detailed comparison to models for a later time (since in de- 
tail, an evolving IMF would have secondary implications for galaxy 
evolution in these models that are not accounted for by this simple 
rescaling). For now, we simply illustrate how galaxy SFRs could be 
significantly lowered by an evolving IMF during the peak epoch of 
cosmic star formation. 

An important constraint on SFR-M* evolution is that its 
amplitude appears t o be essentially con stant from 2 ~ 2 — 
6 ( Stark etal. 2002; ISouwens etdl 12011) . Simulations, on the 
other hand, tend to predict an amplitude that continues to rise, al- 
be it at a slowe r pace than at lower redshifts. Hence as discussed 
in IPava ( 1201 ih . if the IMF varies in some systematic way with 
redshift in order to reconcile models and data, then the evolution 
must reverse at higher redshifts. This would require that the typical 
SFRs of galaxies at z > 2 should be lower at earlier epochs. Such 
an evolution is at least qualitatively consistent with the observation 
that L* decreases going to higher redshifts (i.e. the Universe is "up- 
sizing" at that early epoch, not downsizing). On the other hand, at 
very high redshifts, the CMB temperature rises accordingly, setting 
a larger minimum ISM temperature (and hence a larger minimum 
Jeans mass). In this sense, it is not quite obvious how to reconcile 
the potential discrepancy between the observed normalisation of 
the SFR- A/, relation and modeled one. A quantitative comparison 
will require a more thorough analysis of the data versus the models, 
which we leave for the future. 



4.4 Sub-millimetre galaxies 

Sub-millimetre galaxies at high redshifts are the most rapidly star- 
forming galaxies known, and hence likely represent an extreme of 
ISM physical conditions where our IMF scenario would predict the 
largest variations. Indeed, theoretical models attempting to repro- 
duce SMGs have sometimes resorted to invoking IMF variations in 
order to explain these systems. 

Semi-analytic models with a canonical IMF have diffi- 
culty matching the number c ounts of high-z in f rared and Sub- 
millimetre luminous gal axies JBaugh et al.l |2005|; lOonzalez et al.l 
l20Ilt iNiemi et al.ll20I2n . and have invoked as a solution that the 
IMF is more top-heavy in these merger-driven starbursting sys- 



tems. On the other hand, the required IMF variation in this sce- 
nario is rather extreme, namely a flat IMF above IMq, which is 
disfavoured by observational constraints on the IMF from measur- 
ing th e gas, stellar, and dynamical mass of SMGs (Tacconi et al.l 
l2008l) . 

The simulations of IPave et alj ( 120101) suggest a different ori- 
gin for SMGs, as the high-mass, high-SFR end of the z ~ 2 — 3 
galaxy main sequence. While these models can broadly reproduce 
many properties of SMGs, they do not achieve the high SFRs in- 
ferred from the SMG observations, falling short by a factor of ~ 3. 
This is similar to the correction factor predicted in Equation|4]for 
the most IR luminous galaxies. Hence if the IMF varied in the man- 
ner considered here, this scenario for the origin of SMGs would be 
strengthened. 

Conversely, recent hydrodynamic simulations of isolated 
galaxies and galaxy mergers (similar to those presented here), cou- 
pled with dust radiative transfer calculations and semi-empirical 
cosmological simulations have shown that the number counts of 
high-z sub-millimetre sele cted galaxies may be acc o unted for while 
utilising a Kroupa IMF (Naravanan et alj 1201 (J : iHavward et al.l 
l20ICtl201ir) . at least within the current uncertainties. These models 
suggest that pairs of galaxies blended within a single far-IR beam 
may be an important observational bias at lower flux densities, and 
that the relationship between sub-millimetre flux density and SFR 
may not be linear owing to substantial variations in dust tempera- 
tures during the merger phase. If the observed number counts from 
upcoming surveys turn out to be in good agreement with these pre- 
dictions at a range of far-IR wavelengths, then an IMF-based solu- 
tion may be unnecessary. 

While they diffe r in detail, both the iDave et al.l ( |2010|) and 
iHavward et al.l ( 120111) models suggest that the lowest luminos- 
ity SMGs (Fg5o ~ 5 mjy) galaxies may represent individual, 
heavily star-forming discs, whereas more luminous SMGs likely 
owe their origin to galaxy mergers. In hydrodynamic cosmo- 
logical simulations, individu al discs rarely sustain SFRs at rates 
larger than r^ 500 M ^yr"^ dPekel et alJl2009bl ; lDave et alj|20ld : 
IHavward et al.ll201 11) . Utilising Equation|4l this corresponds to lu- 
minosities Lboi ~ 10^'' Lq , and may reflect a transition lumi- 
nosity where objects above this are typically mergers at z ~ 2. 
Recent analysis of the specific star formation rat e s and Lir- 
Tdust relation for high-z SMGs by iMagnelh eT^ ( |20I2[) sug- 
gests that this may indeed represent a transition luminosity be- 
tween galaxy populations. At Lboi ~ 10^^ Lq , the LiR-Tduat 
relation steepens dramatically, and galaxies rise abov e the "main- 

2 dPaddi et al. 



seque nce" of specific star fo rmation rates at z 
2007i: lRodighiero et akluOIlf) . Moreover, models bv lHopkins et al. 



(2010y show that Lboi ~ IO^^Lq represents a transition lumi 
nosity above which the luminosity function is dominated by merg- 
ers at z - 2. If a Salpeter IMF holds f or z - 2 SMGs, a SFR of 
500 M0yr~^ would correspond to '^ 3 x IO^'^Lq , an implausibly 
low transition luminosity given z ~ 2 merger ra tes (IFakhouri & Mai 



l2008l; IGuo & Whitel2008l:lHoDkins et alj|2010l) 

In short, achieving the extreme SFRs inferred for SMGs at 
their observed number densities using a canonical IMF remains 
a challenge for hierarchical-based models. While observations are 
sufficiently uncertain to be broadly consistent with a variety of sce- 
narios, tensions could generally be alleviated if the SFRs were a 
factor of several lower in these systems as would be the case for the 
varying IMF presented here. Our model form for the IMF suggests 
a reasonable transition luminosity for the break between discs and 
mergers at z '^ 2, given a rough maximal SFR for a quiescent disc 
of ~ 500 M0yr-i (Lboi ~ lO'^L© ). 
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Figure 6. Evolution of c osmic star form a tion ra te density. Purple solid line 
shows model form from iHopkins et al.l yOlOj), assuming a Kroupa IMF. 
This provides a good fit to the observed SFR density measurements of 
(Hopkins & Beacom 2006). When applying our model form for IMF vari- 
ations, the SFR density decreases, as is shown by the red dashed line. The 
blue dash-dot line shows the Wilkins et al. (2008) SFR density that would 
be necessary to match the observed present-day stellar mass density. The 
application of our model IMF decreases the inferred SFR density. But be- 
cause much of the star formation occurs in lower-luminosity systems where 
the IMF is not substantially different from canonical, the amount of the 
decrease is not quite enough to account for the discrepancy between the in- 
tegrated SFR density and stellar mass density as presented in Wilkins et alj 
120081) . 



4.5 The Evolution of Cosmic Star Formation Rate Density 

Deep surveys have permitted measurements of the evolution of the 
star formation rates per unit volume over cosmic time. Similarly, 
increased sensitivities in the near-infrared have allowed for a more 
robust reconstruction of the global stellar mass evolution in galax- 
ies out to 2 ~ 3. In principle, the integral of the star formation rate 
density, corrected for stellar mass loss, should be equivalent to the 
stellar mass density at a given redshift. A number of authors have 
noted, however, that there may be a discrepancy between these two 
quantities, such that the cosmic SFR density appears to be over- 
estimated (Hopkins & Beacoml l2006l : |Per ez-Gonza lez et al.ll2008l : 
IWilkins et al. 2008) . As alluded to in §[T1 an IMF weighted toward 
more high mass stars at high- 2: may go some length toward solving 
this issue. 

Does an IMF that varies with the Jeans mass of GMCs re- 
duce the inferred SFRs of high- 2: galaxies enough to remove the 
discrepancy between the integrated SFR density and the observed 
stellar mass density? To answer this, we need to determine how 
a varying IMF would impact the inferred global SFR density. To 
do this, we employ a functional form for the observed luminosity 
function of galaxies as a function o f redshift that is quite success- 
ful in matching existing data from .Hopkins et al.l ( 12010.) . We inte- 
grate this functional form as a function of luminosity to calculate 
the typical number of galaxies per Mpc^ as a function of redshift 
(integrated to O.IL*). We similarly utilise this information to cal- 
culate a luminosity and SFR density as a function of redshift. This 
information combined with the average number density of galax- 
ies as a function of redshift allows us to calculate the typical SFR 
per galaxy as a function of redshift. We then apply the correction 



giv en in Equation [4] a fter converting to an infrared luminosity via 
the lKennicuta (1998a ) relations, and assuming Lir ^ Lboi- 

Figure |6] shows the resulting SFRD evolution assuming 
an unevolving Kroupa IMF (purple), versus the model evolv- 
ing IMF assuming Mc oc Mj (red dashed). The unevolving 
IMF curve is in goo d agreement with the SFRD compilation by 
iHopkins & BeacomI (2006) when corrected to a Kroupa IMF; for 
clarity, we do not show the individual data points. While the evolv- 
ing IMF makes a large difference to the SFRs of very rapidly 
starforming galaxies, the global effect is less dramatic, since the 
SFRD tends to be dominated by more modestly starforming sys- 
tems. Nonetheless, it is clear that it lowers the SFRD significantly 
at all 2 > 1^ 

.Wilkins et al.l ( 120081) suggested a variable IMF that was em- 
pirically constructed to reconcile the cosmic SFRD evolution with 
the present-day stellar mass density. The resulting SFRD evolu- 
tion in this model is shown as the blue dot-dashed line in Figure|6] 
The IMF model presented here results in less of a correction to 
the present-day stellar mass, since it is clearly closer to the origi- 
nal unevolving case over much of cosmic time. Hence it remains to 
be seen if the IMF model in this work can actually reconcile these 
measurements, although it is clear that it goes some distance to- 
wards this. The standard Kroupa IMF forms roughly twice as many 
stars today as are accounted for in stellar mass density measure- 
ments, whereas our model IMF forms 1.3x the Wilkins et al. esti- 
mate. 

This said, an evolving IMF is not necessary to reconcile the 
differences between the cosmic SFRD evolution and stellar mass 
density. For example, utilising cosmological SPH simulations, ? 
suggest that the two can come into better agreement if the stellar 
mass density measurements are missing faint, low-mass galaxies. 
Similarly, ? have suggested that the evolution of the stellar mass 
density and integral of the cosmic star formation rate density may 
not be discrepant. They argue that if one integrates the luminosity 
function to a common limit (accounting for the evolution of L* 
with z), then some of the observed discrepancy between the star 
SFRD and stellar mass density can be removed. 

Overall, the global SFR density is less affected by our model 
IMF than the most rapidly star-forming galaxies, since much of the 
global star formation occurs in modestly star-forming systems. Our 
model IMF goes some of the way towards alleviating the putative 
discrepancy between the integrated cosmic SFRD and the global 
stellar mass evolution, although current observations still disagree 
on the magnitude of the discrepancy. 

Finally, we note that the rising CMB temperature with redshift 
sets an increasing minimum Jeans mass that is non-negligible at 
large (z >3) redshifts. However, due to the rising SFR per galaxy 
with redshift, the average conditions are already warm enough due 
to gas-dust coupling that the effect of the CMB is negligible. As 
an example, at 2 ~ 4, the CMB temperature is ~ 15 K. The mean 
SFR per galaxy is ~ 20 MQyr^^ , which corresponds to a typical 
gas temperature of '--^ 20 K (see e.g. Figure [T]). 



5 DISCUSSION 

5.1 Arguments for and against a varying IMF 

Probably the single greatest piece of evidence for an invariant IMF 
is that no conclusive evidence for a varying IMF has ever been seen, 
even when probing a fairly wide range of physical conditions and 
environments in star forming regions within the Milky Way and 
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nearby galaxies. For this reason as well as others, all arguments for 
a varying IMF must b e taken with a healthy degree of skepticism. 
The recent review by ISastian et al.l (120100 takes a rather critical 
view of observational claims of IMF variations; we refer the reader 
there for a fuller discussion. 

There are many theoretical studies that suggest th at the IMF 
should be essentially constant in most relevant regimes. iKrumhola 
( 12011 ) present a model in which the temperature in typical star 
forming regions in GMCs is set by radiation feedback from ac- 
cretion onto protostars. In this scenario, one can aiTive at a charac- 
teristic stellar mass primarily via fundamental constants, with only 
weak scaling with the physical properties of the GMC. Similar ar- 
guments were made by Bate! ( l2009h . This theory may explain the 
near constancy of the IMF in normal star forming regions in the 
Galaxy, and has the additional attractive aspect of explaining why 
stars form to massive enough scales to begin hydrogen fusion. 

On the other hand, there have been a number of recent claims 
that the IMF varies with galaxy properties. It is now reason- 
ably well established that the H q to UV luminosity ratio varie s 
syste matically with luminosity JHoversten & Glazebrookl l2008h. 
SFR JLee et alj 1200 9) and/or surface brightness (IMeurer et al.l 
l2009l) . in the sense that for fainter galaxies there is less Ha than 
expected given the UV luminosity. Since Ha traces very high mass 
stars while UV traces moderately high-mass stars, one possible ex- 
planation is that smaller GMCs are unable to produce very high- 
mass stars. In this scenario, fainter systems with on average smaller 
GMCs would have an integrated galactic IMF that is less top-heavy. 
This variation would work in the same direction as the IMF vari- 
ations considered in this paper (i.e. more top-heavy/bottom-light 
in more rapidly star-forming systems), but for a different reason; 
also, it is seen to operate down to very low SFRs (and indeed it 
is most evident there), whereas our results suggest that IMF varia- 
tions based on the Jeans mass conjecture have little impact at SFRs 
of the Milky Way or smaller That said, numerous investigations 
recently have discussed other physical effects that could lead to the 
observed trends in HaAJV, such as the leakage of ionizing photons, 
stochasticity in the formation of the most massive stars (?), variable 
star formation histories (?), and extinction effects. Similarly, stud- 
ies investigating the production rate of ionising photons in stellar 
clusters (?), as well as outer galaxy discs (?) find no evidence for 
a varying IMF. Hence while intriguing, Ha/UV variations are not 
widely considered to be an unambiguous indication of a varying 
IMF. 

E arly-type galaxies also show ta ntalizing hints for IMF vari- 
ations. |vm^o^uin&_Conroi] (|201Jj) observed the Wing-Ford ab- 
sorption bands (tracing stellar masses < 0.3Mq) in nearby cluster 
elliptical galaxies, and found that they indicate a much larger popu- 
lation of low-mass stars than expected from a Galactic IMF, i.e. tha t 
the IMF in these galaxies is bottom-heavy. ICappellari et alj ( l2012n 
did careful dynamical modeling of ellipticals from the ATLAS-3D 
sample, and showed that high-mass ellipticals tend to have such 
mass-to-light (M/L) ratios that they can only be explained by a 
bottom-heavy IMF. They also found a strong trend of IMF variation 
with M/L. These data argue for a bottom-heavy IMF as opposed to 
a bottom-light one as considered in this paper, but a naive extrapo- 
lation of their trend with M/L suggests that rapidly star-forming 
systems may have a top-heavy/bottom-light IMF. However, it is 
also expected that massive star- forming systems at high-z should 
evolve into present-day passive ellipticals, and thus it is not imme- 
diately clear how to reconcile these differing IMFs within a single 
galaxy evolution scenario. If, indeed, the IMF is bottom-heavy in 
the progenitor galaxies of present-day ellipticals, then this likely 



rules out the model p resented in this paper (Ivan DokkumI BOOSt 
ICappellari et"Zll2012l; ?). 

If the IMF is indeed invariant, then what of the cosmic star 
formation problems discussed throughout this work? There are 
certainly solutions that do not require an IMF-based solution. 
There are many well-documented systematic uncertainties associ- 
ated with measuring SFRs and Af,, particul arly at high redshift s 
when only broad-band data is available (e.g. IShaplev et al. 2005]). 
It is not difficult to imagine that one or more of t hese could a lleviate 
the apparent discrepancies (see a discussion in iDavauOOSO . There 
could also be failures in current modeling techniques or assump- 
tions. For instance, the high sSFR of z ~ 2 galaxies could owe 
to a strong inefficiency of star formation in earlier galaxies owing 
e.g. to lower metallicities JGnedin & Kravtsovll201Cll ; lBolatto et al.l 
uOllf) . resulting in an acc umulated gas reservoir tha t fuels excess 
star formation at z ~ 2 dKrumholz & DekelluOllh . This would 
not alleviate differences between the observed cosmic SFR and 
h'h buildup, but perhaps those could be explained by other sys- 
tematics; for instance, if the bottom-heavy IMF seen in early-type 
galaxies is real, the total cosmic stellar mass in this component 
may be underestimated. Finally, it is worth noting that numeri- 
cal unce rtainties still abound in modeling galaxy formation: for 
instance, iKeres et alj ( 120111) found that SFRs in (~ L* ) galaxies 
modeled with the moving-mesh code AREPO are '^ 2 — 3 greater 
than those modeled using the SPH code GADGET-3. With all the 
various factor-of-two uncertainties in play, it is difficult to argue 
strongly for IMF variations as being the only viable solution to any 
given problem. 

As discussed in § 14.11 a varying IMF may bring the ob- 
served Kennicutt-Schmidt relation index down from ~ 2 to ~ 1.5, 
in agreement with most models which suggest a star formation 
rate regulated by dynam ical processes. However, as shown by 
lOstriker & Shettvl (1201 ih , if the star formation rate is instead regu- 
lated by supemovae-driven turbulence in a medium where the verti- 
cal pressure is dominated by gas, then one would expect a KS index 
of ~ 2. In this sense, even the reduction of the observed KS index 
is not reason enough to require a variable IMF. 



5.2 Arguments for and against the Jeans mass conjecture 

In this work, we make the assumption that the IMF characteristic 
mass scales as the GMC Jeans mass. Here we discuss arguments in 
favor of and against this Jeans mass conjecture. 

Stars form in dense cores in GMCs (e.g. Evansl 1 19991 ; 
iLada & Ladal 120031) that ap pear to have a ma s s function similar 
in shape to the stellar IMF (lAndre etal.ll 19961; lMotteetaLll200ll; 
I Johnstone et"ai]|200ll ; lNutter & Ward-Thompsonll2007h . This sug- 
gests that stars acquire their mass from the cores they are born 
in, and are thus plausibly dependent on the fragmentation scale 
of the cloud. In the same vein denser(less-dense) clouds, which 
would have a lower(higher) Jeans mass given the Mj ~ T^'"^ /ii}''^ 
scaling, should have a more bottom-heavy(bottom-light) IMF. Ob- 
servations of the Taurus star-forming region and Orion Nebula 
Cluster provide some evidence of this effect (IBriceiio et alj|2002l ; 
iLuhman et al.ll2003l ; lLuhm"ai]|2004l ; ?). 

Numerical calculations also suggest that the number of frag- 
ments that form in a GMC is related to the number of Jeans masses 
available tBate & Bonnell 2005). Klessen et al. (2007) showed via 
numerical simulation that warm gas comparable to what may be 
found in nuclear starbursts indeed forms a mass spectrum of col- 
lapsed objects with a Mc~ 15 Mq, though the degree of frag- 
mentation can depend on the exact form of the equation of state of 
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the gas. It should be noted that observations of the Arches (?) and 
Westerl und 1 clusters (??), r egions similar to those in the simula- 
tions of iKlessen et alj 120071) show minimal evidence for a varying 
IMF.^ 

lElmegreen et alJ ( 120081) asked, given cooling rates by molec- 
ular lines, and heating rates by dust-gas energy exchange in dense 
gas, how should A/,j scale with density? These authors find a weak 
dependence of the Jeans mass on the density as n^'*. Examining 
Figure[T] the mass-weighted mean density of GMCs is ^ 10* cm^'^ 
during the starburst, compared to a few xlO^cm^'^ during quies- 
cent mode. Following the Elmegreen et al. scaling^ this alone 
would lead to an increase in Jeans mass of ~ 4. In other words, 
even with a relatively weak scaling of Mj with density, the condi- 
tions inside nuclear starbursts are extreme enough to result in vari- 
ations in the GMC physical properties as found in our models. 

A more general argument against the conjecture that A/c oc 
M} is that the choice for the density that sets the Jeans mass is 
not obvious. In this work we sidestep this problem by consider- 
ing only how Mj scales with galaxy properties and assuming that 
Mc oc Mj without worrying about the constant of proportional- 
ity. However, using typical numbers from our simulations results 
in Mj ~ 25OAf0 for the Milky Way (Figure [2) which is un- 
satisfyingly far from Afc ~ O.SAffn. lLarsonl OOOSD assumes that 
the characteristic density is set by when the gas becomes thermally 
coupled to the dust, and argued ag ainst this scale be ing set by tur- 
bulent or magnetic pressure. But in Krumhola ( 1201 ih and other re- 
cent models, the density where isothermality is broken is set in- 
stead by radiative feedback, and thus the IMF characteristic mass 
is not directly related to the Jeans mass. We note that even if the 
A/c is set by radiative feedback, the interstellar radiation proper- 
ties vary substantially among our various simulated galaxies, and 
hence there may still be some IMF variations predicted in such a 
model. Nonetheless, the physical origin of the IMF characteristic 
mass remains uncertain, and it is by no means entirely clear that it 
is directly related to the typical GMC Jeans mass. 

In this paper we suggest that heating of H2 gas by young 
stars formed during rapid star formation may increase the Jeans 
mass. Other mo dels find simila r result s, though via different heat- 
ing mechanisms. IPapadopoulod ( 12010 ) find that the cosmic ray en- 
ergy densities in starburst environments may be increased, lead- 
ing to warmer gas kinetic temperatures of order ~ 80 — 160 K. 
IPapadoDOuIos et alj ( 1201 In find that this can lead to increased Jeans 
masses up to a factor ~ 10. The relevant physical processes in 
the Papadopoulos model are different from those presented in our 
work. In a scenario where cosmic rays dominate the heating, the 
gas and dust do not necessarily have to be thermally coupled. Nev- 
er theless, the results are qui te similar to what is found in our study. 
In [Narayanan et al.l ( 1201 id) , we found that even when assuming a 
cosmic ray flux that scales linearly with the star formation rate, the 
heating of dust (and subsequent dust-gas energy exchange at high 
densities) by young stars in starbursts tends to do minate the heating 
of the gas. Similarlv. lHocuk & Spaana ( 120100 and lHocuk & Spaanj 
( 120 llh have suggested that a strong X-ray field in the vicinity of 
an accreting black hole may impact the fragmentation length of 
GMCs, weighting the clump mass spectrum toward higher masses. 
This scenario is likely to be most relevant in rather extreme cases 



"^ Note that a similar scaling relation comparing to the global SFR is non- 
trivial to derive. This depends on radiative transfer processes as well as 
potential contribution to the radiation field from old stellar populations. 



(e.g. near an AGN), and not necessarily representative of quies- 
cently star-forming galaxies. 

Overall, whether the stellar IMF varies in space and/or time, 
and how remains an open question. Strong IMF variations are 
clearly disallowed by observations, but mild variations may still be 
accommodated within cuiTent constraints. The Jeans mass conjec- 
ture is by no means unassailable, and the resulting mild variations 
in the IMF are probably difficult to exclude given cuiTent data. 



6 SUMMARY 

We use hydrodynamic galaxy evolution simulations and sophisti- 
cated radiative transfer models to study how the IMF might vary 
with galaxy properties, under the assumption that the IMF charac- 
teristic mass scales with the H2 mass-weighted Jeans mass in the 
ISM. We find that such an assumption results in an typical Jeans 
mass that varies mildly with star formation rate, roughly as A/j oc 
SFR"'* at star formation rates above that of the Milky Way. At 
SFRs comparable to and lower than that of the Milky Way, there 
is little variation in Af j in our model because cosmic ray heating 
begins to dominate over dust heating, setting a temperature floor 
at ~ 8 K. While this variation is not strong, it is still substantial 
enough to impact the inferred properties of rapidly star-forming 
galaxies. 

We parameterise the implied IMF variation by considering its 
impact on SFRs derived from bolometric luminosities. We find that 
the relation is well-fit by SFR= (Lboi/IO^^L©)"****, as opposed 
to a linear relationship between these quantities as for an invari- 
ant IMF (assuming that any AGN contribution has been excluded). 
Since all extragalactic SFR tracers attempt to measure Lboi, either 
directly in the infrared, via a dust-corrected UV continuum, or via 
nebular emission line measures, this would imply that their inferred 
SFRs should be lowered in more luminous (i.e. more rapidly star- 
forming) systems. For instance, in sub-millimetre galaxies that are 
the most bolometrically luminous systems in the Universe, the in- 
ferred SFRs would be up to a factor of five lower than implied by 
canonical conversion factors. Large, rapidly star-forming disks at 
high redshifts would have their inferred SFR lowered by a factor of 
two. 

We study the impact of such an IMF variation on a variety of 
galaxy properties across cosmic time. These include the star for- 
mation law, the SFR— A/, relation (i.e. the galaxy main sequence), 
sub-millimetre galaxies, and the evolution of the cosmic SFR den- 
sity. In each case, previous work has highlighted tantalizing albeit 
inconclusive evidence for discrepancies between models and data, 
or between star formation rate and stellar mass data. The discrep- 
ancies (if real) all point in the same direction, in the sense that the 
SFR infeiTed from high-mass tracers seems to exceed the observed 
growth in stellar mass by '^ x 2 — 3 during the peak epoch of cos- 
mic star formation. In each case, the varying IMF based on the 
Jeans mass conjecture tends to go towards reconciling the possible 
discrepancies. We caution that there is no smoking gun evidence 
for IMF variations, and there are numerous ways these discrepan- 
cies may be alleviated without varying the IMF. Nonetheless, it is 
interesting that the ansatz of M^ oc Mj makes measurably differ- 
ent predictions for the evolution of galaxies across cosmic time in a 
manner that, at face value, tends to reconcile any possible discrep- 
ancies. 

We remain agnostic on the validity of our underlying assump- 
tion, namely the Jeans mass conjecture, and we do not assert that 
IMF variations are the only solution to the problems discussed here. 
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Rather, this paper merely explores the ansatz that Mc oc Mj, and 
then examines its effects on the star formation rates across cosmic 
time. As we discussed in ^ cuiTent models for the origin of the 
IMF characteristic mass are mixed on whether this assumption is 
correct; a recent class of successful models argues that Ale is set 
by radiative feedback in the ISM, not Mj. We leave further ex- 
plorations along these lines to a different class of models that fo- 
cus on the origin of the IMF. Additionally, such an evolving IMF 
would have a wider range of impacts than explored here, owing to 
differences in e.g. stellar mass return, metallicity, and population 
synthesis that have not been fully accounted for here. Also, at a 
constant ISM density, the rising CMB temperature will force in- 
creasingly higher Jeans masses with redshift, further complicating 
matters. Investigating these effects would require running simula- 
tions that self-consistently model all these processes during their 
evolution, which is something we leave for the future. 
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APPENDIX A: A JEANS MASS THAT VARIES WITH THE 
MEDIAN DENSITY 

It is conceivable that one could consider the density term in the 
Jeans mass equation as the median density (i.e.the density above 
which half the mass resides), rather than the mean density. As a re- 
minder, while we resolve the surfaces of GMCs, we do not resolve 
their internal structure. We assume that our GMCs have a lognor- 
mal density distribution. In this case, the median density can be 
expressed in teiTns of the width of the lognormal, ap] 
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where n is the mean density of the GMC. Numerical ex- 
periments suggest that a'^ ~ ln(l -I- 3M^/4), wh e re M 
is the on e-dimensional Mach number JOstriker et akl uOOlt 
IPadoan & Nordlund 2002 1'^l In Figure lAll we show the relation- 
ship between SFR and Mj when the Jeans mass is calculated using 
Wnicd instead of n. 

The relation between Mj and SFR in Figure IaTI is relatively 
tight. When comparing to Figure|2] we see that the scatter is greatly 
reduced when considering a Jeans mass calculated with the rimed, 
rather than n. In Figure |2] the bulk of the scatter arises from rel- 
atively active systems. For example, at a SFR of 50 MQyr^^, the 
galaxies with the highest average Jeans mass in Figure |2] are on- 
going mergers with warm ISMs, whereas the galaxies with lower 
Jeans masses are either mergers during a more quiescent phase, or 
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Figure Al. Relationship between SFR and Jeans mass when A/j is calcu- 
lated using the median density in GMCs, rather than the mean. The relation 
is not dissimilar from using the mean density. See text for details. 



discs. Using the median density rather always results in a Jeans 
mass that is lower than when using the mean density, due to the 
extra term that scales with M^. Because mergers and heavily 
star-forming discs tend to have high- velocity dispersion gas (e.g. 
[Narayanan et alj|201 Ibl) . their Jeans masses are reduced substan- 
tially when using rimad- This tightens the dispersion dramatically 
in the range SFR=1Q - 100 M0yr"\ 

One consequence of using nmed rather than n is that a signif- 
icant increase in the characteristic mass of the IMF is only seen in 
heavily star-forming galaxies (SFR > lOOMQyr"^). This suggests 
that such a variation may only be seen during short-lived phases 
in local mergers, or in more heavily star-forming galaxies at high 
redshift. 

The results of this paper change little when utilising rimed 
rather than n in the Jeans mass calculation. The normalisation of 
the Mj-SFR relation in Figure IaTI decreases, while the exponent 
increases slightly. These effects cause the normalisation of Equa- 
tion|4]to decrease by a factor of ~ 2, while the exponent decreases 
to ~ 0.85. 



* We note that other a uthors have found a range of possible forms for the 
dependence of o-p. See lLemaster & Stonel 420081) and lPrice et al.l 1 1201 ID . 
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